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^ ; 1 INTRODUCTION 



ABSTRACT 

We create new U, B, V, R and I-band light curve templates of type la Supernovae 
(SNe la) and re-analyze 122 nearby (redshift < 0.11) SNe l a using a new "Multi- 
band Stretch method," which is a revised Stretch method (cf. [Perlmutter et al.lll997t 
iGoldhaber et al.|[200l[ ) extended to five bands. We find (i) our I-band template can 
fit about 90% of SNe la I-band light curves, (ii) relationships between luminosity, 
colours and stretch factors, (iii) possible sub-groups of SNe la, and (iv) the ratio of 
total to selective extinction R in other galaxies can be consistent with that in the 
Milky Way under the assumption that SNe la have diversity in their intrinsic colour. 
Based on these results, we discuss how to select subsets of SNe la to serve as good 
distance indicators for cosmology. We find two possibilities: one is to choose "BV 
bluest" (—0.14 < (B — V) max ^ —0.10) objects and the other is to use only SNe la 
which occur in E or SO galaxies. Within these subsets, we find the root mean square 
(r.m.s.) of peak B-band magnitudes is 0.17 mag ("BV bluest" sample) and 0.12 mag 
(E or SO sample). 

Key words: supernova: general. - light curve. - distance scale. 



SNe la are some of the most important objects in observa- 
tional cosmology because of their role as distance indicators. 
SNe la are bright enough to be observ ed at high redshifts , 
they have uniform peak luminosity (cf lBranch et al.l 1 1992T ) . 
and they are expected to evolve less fr om low z to high z 
than other objects such as galaxies (cf. iRiess et al.ll"l999bh . 
Many studies have used SN la as distance indica tors to show 



and Tenagra Observatory Supernova Search es (LOTO SSQ]), 
the N earby Supernova Factory (SNfactory, lAldering et al.l 
[200 ^, the SDSS-II Sup e rnova Survey (S PSS SN Sur- 
vey, ]Frieman etHI |200S| ; ISako et al.l l2008t ). Equation of 
State: SupErNovae trac e Cosmic Expansion (ESSENCE, 
IWood-Vasev et aL 20071). and th e SuperNova Legacy Sur- 



the accelerating expansion of the universe (e.g Riess et al 
199S : iPerlmutter et al.lll999l; iKnop et al.l 120031: biess et al 
2004 : lAstier et all 120061 : IWood-Vasev et all 120071 ). All of 



Ham uv et al 
iReindl et al 



these works are based on the empirical relati onship betwee n 
light curve shapes and peak luminosity (cf. |Phillipslll993l ). 
From the 1990s onward, as more low-redshift SNe la have 
been found, the relation has been revised (cf, 
19961 : IRiess et al.lll999al : lAltavilla et~al1 12004 
2005l h 

Recently, the number of well-measured SNe la has 
grown larger and larger, especially at intermediate and 
high redshifts. There are many supernova survey pro- 
grams currently running; for example, the Lick Observatory 



vey (SNLS, iPritchet et all l2005r i. Most of SNe la found 
in these surveys are observed in multiple passbands and 
over many epochs. These high quality samples have given 
us clue s to underst a nd the nature of SN la. For ex- 
ample, lHamuv et alj (fl99 6b) showed that brighter SNe 
la, like SN1991T, tend to occur in galaxies with active 
star formation, while fainter SNe la, such as SN1991bg, 
ten d to occur in galaxies which lack star formation (see 
also [ Howell 200 l|; Ivan den Bergh et al.ll2005l : ISullivan et ail 
2006). Benetti et al.l l|2005l ) discussed the intrinsic diversity 
of the spectroscopically normal SN la (called " Branch nor- 
mal", ^^^^^\^9^).^^^^^^^^i) pointed out 
there are two subgroups within Branch normal SN la. Un- 
usual SNe la have also been found: SN2000cx and SN2002cx 
were reported as an extreme samp le of spectroscopically pe- 
culiar SN la (|Li et al.ll200ll . l2003l h; iHowell et all <|2006h re- 
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Table 1. Number of SNe la in each band 



Sample 


U-band 


B-band 


V-band 


R-band 


I- band 


1A 


49 


122 


122 


101 


115 


IB 


41 


108 


108 


88 


102 



1A are all SNe la we analyzed in this paper. 
IB are selected SNe la based on x 2 ( see §3.2). 



ported SNLS-03D3bb was a Super-Chandrasekhar-mass su- 
pernova with an exceptionally high luminosity and low ki- 
netic energy. Ha emission in the spectra of SN2002ic and 
SN2005gj at early phases indicates the e xistence of circum- 
steller material aro und the progenitor l|Deng et al. 1 12004 
lAldering et alj|2006l and so on). 

Although SNe la seem to be good standard candles for 
cosmological studies, there are several issues we must under- 
stand in order to use them optimally. The main issues are 
(1) intrinsic photometric and spectroscopic diversity of SNe 
la and (2) extinction in host galaxies. Both issues severely 
limit the use of SNe la as a standard candles if we can not 
deal with them properly. Recent studies h ave revealed an 
intrin sic di versity of SNe la; fo r example, iMannucci et al.l 
(|2005h and I Sullivan et all (120061 ) showed there are two types 
of SN la, called "prompt" and "delayed " , which differ in the 
distribution of their stretch factors, and lBerian et al] (|200d ) 
demonstrated how spectroscopic diversity can affect signif- 
icantly the intrinsic colours of SNe la and K-corrections. 
The latter issue is important if we want to use larger SN la 
samples. Correction of extinction by dust in the host galaxy 
is necessary because almost all SNe la may be affected by 
the dust. Previous works insisted that the ratio of total to 
selective extinction, Rb, i n other galaxie s is sm aller than 
that of t he Milky Way: see lAltavilla et alj (l2004h (her eafter 
ALT04) , iPhillips et all (|l999l ). and lKnop et all feoOSl ) . 

In this paper, we discuss these issues using a sample of 
nearby SNe la. We re-analyze published data on nearby SNe 
la with an improved "Multi-band Stretc h method" which 
is based on the usual "stre tch method" l|Perlmutter et all 
ll997l ; lGoldhaber et al.|[200ll ). 



2 PHOTOMETRIC DATA 

We use U, B, V, R and I-band light curves in Vega sys- 
tem of 122 SNe in this work. All of the sample have 
B and V-band photometry, but they don't always have 
U, R and I-band photometry (see Table [TJ. All photo- 
metri c data are obtain ed from published pap ers, mainly 



from 



l|l999al ) and 



Jha et all (l2006h ( h ereaft er 
Hamuv et al.1 Jl99fj). 



JHA06), iRiess et al.l 
We also took photo- 



metric data fro m iFord et all (119931), I Wells et all |jl994l) 
iRichmond et all (I 19951 ), iPatat et all ill99rj), | Lira et al.l 



(1998), 



Krisciunas ct al. 

17 



(2000), 



iSuntzeff et all J 19991), iKrii 
iModiaz et al.1 (l200ll). 'iHo et al.1 120011), IKrisciunas et al 
<l200ll). IStrolger et al.l d2002f). IKrisciunas et al.1 Em, 
Li et all ll2003).IVinko et al.l d2003l) IValentini et al.l 1200. 



.19991), 
(l200ll). iHoet al.1 



Vinko et al.l ll2003l),lGarnavich et all (|2004l).lGermanv etall 
<|2004|). IKrisciunas et all ll2004al). IKrisciunas et~all (l2004bl), 
iPignata et al.l (|2004l ), lAnupama et al.l l|2005l ). iKuntal et all 



l|2005l ) and lSahu et~ai1 l|2006l ). All of SNe la used in this pa- 
per have B and V-band data. We analyzed only bands in 
which there are more than 3 measurements. For example, if 
a SN had five B-band measurements, five V-band measure- 
ments and two I-band measurements, we did not include 
its I-band photometry in our analysis. There are some SNe 
la which were observed by different groups independently, 
but we didn't combine these light curves to avoid systematic 
errors; instead, we used the light curve with the largest num- 
ber of epoc hs. Many of our data are al so found in the com - 
pilations of lAltavilla et al.1 (|2004l ) and iReindl et al.1 (|2005l ). 
Figure [T] shows the redshift distri b ution of our sample. 

According to IReindl et al.1 (|2005l ), there are five 
"SN1991T-like" SNe (SN1991T, SN1995ac, SN1995bd, 
SN1997br, SN2000cx) and seven "SN1991bg-like" SNe 
(SN1991bg, SN1992K, SN1997cn, SN1998bp, SN1998de, 
SN1999by, SN1999da) in our sample. The other 111 SNe 
are either "Branch normal" or lack detailed spectroscopic 
analysis. 



3 ANALYSIS 

3.1 Light Curve Parameterization 

There are several ways to describe the shape of a super- 
nova's light curve, but we can classify them roughly into 
two types. One type concentrat es on the dat a in a single 
passband, such as the " A ml5" (|Phillipslll993l ) or "Stretch" 
l|Perlmutter et alii 19971 ) algorithms. The other type includes 
measurements of light curves in multiple passbands, such 
as th e "MLCS, M LCS2k2" (IRiess et al.l Il99rj ; |jha et al.1 
120071), "CM AGIC " J Wang et al.ll2003l )."arid "SALT, SALT2" 
( Guv et al.ll2005l , l2007l ) algorithms. The latter type is more 
complex and empl oys SN colours as a part of the analysis. 
IPrieto et alj (|2006h proposed a new method, a combination 
of A ml5 and MLCS, to gain the advantages from both 
approaches. 

For this work, we chose the single-band stretch method 
as a base, since it is simple to parameterize light curves: 
we just need one template for each band. We created new 
U, B, V, R and I-band templates (see §3.4), which we sim- 
ply stretch to fit observed light curves, without having to 
account for additional assumptions such as dust extinction. 
This method works well when we simply want to assign a 
light curve shape and peak apparent brightness in each band. 

Usually, the stretch method has been applied only to 
B or V-band light curves. However, in this work, we apply 
the stretch method to data in all passsbands: U, B, V, R 
and I. In other words, we fit up to 11 parameters to each 
event: U, B, V, R and I-band stretch factors, U, B, V, R and 
I-band peak magnitudes, and a time tsmax for maximum 
light in the B-band. The reason that we apply a stretch 
factor to U, B, V, R and I-band in dependently inst ead of a 
common stretch factor like SALT ([Guv et al.N2005l ) is that 
we wish to investigate the diversity of light curve shapes 
between passbands. We call this method the "Multi-band 
Stretch method" . 

It is well known that there are some peculiarities in 
redder band light curve shapes, especially in I-band. There 
are usually two peaks in I-band light curves, but some SNe 
la lack the second peak. There are very few of these SNe 
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la without the second peak in our sample, so we use only 
"normal" I-band template for our analysis (see §4.3). 

Using information in multiple passbands improves the 
determination of the time of B-band maximum brightness. 
By using data from the second peak in R and I-band around 
30-40 days after B-band maximum, we can determine the 
time of B-band maximum even if there are no observations 
around the first peak. In addition, as a result of fitting, we 
can use colour information to study the intrinsic colours of 
SNe la, the effects of dust extinction, and so on. 

Please note that we use the inverse of the usual stretch 
factor (1/s) for easier comparison with the Aml5 parameter. 



3.2 Fitting Algorithm 

We used a reduced \ 2 method for fitting templates to ob- 
servations. Our definition of reduced x~ ls as follows; 



2 

X = 



DOF 




(Am : 



mxi 



K 



X,(ti 



»)/(!+*) 



M^ mp {(l+Z)x(l/S (x) ) X (ti-tBmax)} 



Xn 



-/<) 



Here, n is the total number of observed epochs, DOF is 
the number of degrees of freedom in the fitting, mxi is the 
observed apparent magnitud e in X-band after corr ection for 
Galactic dust extinction from lSchlegel et all (|l998T l. t Bm ax is 
time of the maximum in B-band, Kx,t is the K-correction in 
the X-banc0 magnitude at time t — tBmax , H is the distance 
modulus derived from a recession velocity, M^ mp (t) is the 
template's magnitude at time t, six) is the stretch factor in 
X-band, X max is the peak magnitude in X-band at the time 
of B-band maximum, a a t s is the photometric error at the 
measurement, and a te m P iate is the error of template at the 
epoch. We fit the parameters tsmax, s^x), X m ax to each set 
of curves of one supernova. Although our Eq.(l) is expressed 
in magnitudes, we performed the fitting in linear scale. U, 
B, V, R and I-band light curves were treated with equal 
weights. We used the recession velocity of each host galaxy 
and cosmological para meters (Hp = 70.8 k m/s/Mpc, Om = 
0.262, fi a = 0.738 from lSpergel et alfeOOTT l to determine the 
distance modulus. Note that we did not introduce any colour 
terms in Eq.(l). We estimated the typical size of errors in our 
fitting procedure with Monte Carlo simulations. We fitted 
30 artificial multi-band light curves which made from light 
curve templates and obtained the dispersion for each fitted 
parameter. 

In our Multi-band Stretch method, the template can 
be applied to the observed light curve anywhere from -10 
days to +80 days after B-band maximum (-10 days to +70 
days for U-band only). These ranges are wider than th ose of 
previous papers. For example, iGoldhaber et al.l l|200lf ) used 
a range of -25 days to +50 days. As a result, we can use 
information from later phases, especially around the second 
peak in R and I-band. 



3.3 K-corrections 

Even though SNe la are simpler objects than galaxies, their 
K-corrections are not so simple since one must account 
for their spectral evolution. In this work, we define a K- 
correction as follows: 

fXF(X,t)Sx{X)dX 
J XF[X/(1 + z),t]Sx(X)dX 

Here F(X, i) is the spectral energy distribution (SED) of 
the SN at time t, Sx (A) is the effective spectral transmission 
in the X-band, and Kx .t is the value of K-correc tion at time 
t. See lKim et ail (|l996l ) and lNugent et all (|2002T l for details. 
We used Bes sell U, B, V, R and I-band filter shapes (c.f. 
iBessell |[l990h for the calculation. 

St udies of K-correction s in nearby SN la samples goes 
back to lHamuv et al.l (ll993l).lKim et al.1 (|l996l ) followed this 
work, and Nugent et al.l ( 20021 ) provided a synthetic spectral 
library for SNe la. We used the Nugent's model versionl.10 
for our calculations of K-corrections. There are three type 
of spectral models, named "Branch normal", "SN1991bg- 
like", and "SN1991T-like" llBranch et al.ll 1993 
classification of each SN in 



K c x °r tS = 2.51og(l + 2 )+2.51og 



(2) 



We use the 



Reindl et al.l (|2005l ) to choose a 



spectral model, or used the "Branch normal" spectral model 
when there was no spectral classification.. 

3.4 Template 

Figure [2] shows our new U, B, V, R and I-band templates of 
SN la light curves. The templates were composed from 110 
nearby (z < 0.05, 90% have z < 0.03) SNe la light curves. 
The recipe to make the templates is as follows: 

(i) Select a de nsely-observed " seed" light curve (we se- 
lected SN2001el (|Krisciunas et alj|2003l) for this work) and 
interpolate epochs. 

(ii) Stretch the other 110 nearby SNe la light curves to 
fit in the current templates without K-corrections, then take 
the weighted average of all. 

(iii) Iterate step (ii) until the templates become smooth. 

(iv) Finally, s tretch each template to fi t in the s = 1 light 
curve defined bv lPerlmutter et al.l l|l997h templatfl 

We didn't apply K-corrections to make the light curve 
template because well-observed SNe la which control the 
shape of light curve template are at the lowest redshift, and 
the size of K-correction is smaller than typical photometric 
errors in all bands. 

We calculated the residuals of each light curve from the 
template (see Figure [2}. We regard this scatter from the 
template as a result of the intrinsic diversity of SNe. The 
templates are available from our web sit^f]. 



2 We didn't need cross filter K-corrections because the typical 
redshift of our sample is so small. 



3 http:/ /supernova. lbl.gov/~nugent/nugent_templates. html 

4 http:/ /supernova. lbl.gov/~nugent/nugent_templates. html 

5 http:/ /www. ioa.s. u-tokyo.ac.jp/~takanashi/works/SN/templatc/ 



Table 2. The Best Fitting Parameters of 108 Selected SNe la 



SN name 


U 


{ U ) 




B 




S(B) 




V 




S(V) 




R 






{H) 




1 








s 


T\ 
1 ) 




x 2 

d.o.f. 


host 


ref. 


SN1989B 


-18.61(0.03) 


0.988(0.018) 


-18 


69 


0.04) 





902 


0.021) 


-18 


98 


0.02) 





999 


;o.oi2) 


-19 


20(0 


02) 





944(0.020) 


-19 


26 





02) 





936 


0.028) 


30.89 


0.476 


- 


(1) 


SN1990O** 


- 


- 


-19 


41 


'0.02) 


1 


109 


0.032) 


-19 


30 


0.02) 


1 


138 


;0.024) 


-19 


29(0 


03) 


1 


104(0.025) 


-18 


95 





02) 


1 


061 


0.026) 


35.60 


0.891 


- 


(2) 


SN1990T 


- 


- 


-18 


92 


0.06) 


1 


067 


0.048) 


-19 


06 


0.06) 





995 


;0.032) 


-18 


95(0 


03) 


1 


050(0.020) 


-18 


74 





03) 


1 


006 


0.017) 


36.18 


0.387 


E 


(2) 


SN1990Y 


- 




-18 


29 


0.07) 


1 


076 


'0.054) 


-18 


67 


0.09) 


1 


Oil 


;0.043) 


-18 


86(0 


09) 





954(0.045) 


-18 


79 





17) 





864 


0.056) 


35.95 


1.445 


E 


(2) 


SN1990af* 


- 


- 


-18 


91 


0.01) 





787 


'0.014) 


-18 


87 


0.01) 





789 


0.011) 




- 






- 




- 












36.69 


0.392 


E 


(2) 


SN1991S 


_ 


_ 


-19 


32 


0.13) 





952 


0.067) 


-19 


02 


0.04) 


1 


096 


;0.034) 


-19 


00(0 


03) 


1 


091(0.021) 


-18 


85 





06) 





996 


0.039) 


36.89 


0.478 


_ 


(2) 


SN1991T 


-21.26(0.02) 


1.155(0.023) 


-20 


77 


0.04) 


1 


119 


0.025) 


-20 


92 


0.02) 


1 


102 


;o.oo2) 


-20 


93(0 


03) 


1 


092(0.030) 


-20 


73 





04) 


1 


103 


0.029) 


32.34 


2.363 


- 


(3) 


SN1991U** 






-19 


26 


0.06) 


1 


005 


0.028) 


-19 


14 


0.06) 


1 


102 


;0.032) 


-19 


17(0 


03) 


1 


117(0.025) 


-19 


12 


;o 


02) 





989 


0.024) 


35.65 


0.427 


- 


(2) 


SN1991ag* 




_ 


-19 


41 


'0.02) 


1 


120 


0.016) 


-19 


38 


0.03) 


1 


110 


;o.oi4) 


-19 


38(0 


03) 


1 


070(0.025) 


-19 


07 


;o 


01) 


1 


065 


0.014) 


33.94 


0.311 




(2) 


SN1992J 


- 


- 


-18 


69 


0.07) 





824 


0.037) 


-18 


89 


0.21) 





784 


'0.111) 




- 






- 


-18 


55 





03) 





804 


0.024) 


36.44 


0.543 


- 


(2) 


SN1992K 


- 


_ 


-18 


39 


0.07) 





590 


0.036) 


-18 


60 


0.23) 





636 


'0.101) 




- 






- 


-18 


17 





04) 





801 


0.022) 


33.37 


0.831 


- 


(2) 


SN1992P** 


- 


- 


-19 


22 


0.02) 


1 


035 


0.020) 


-19 


11 


0.01) 


1 


044 


;o.oi6) 




- 






- 


-18 


75 


;o 


03) 


1 


127 


0.057) 


35.27 


1.358 


- 


(2) 


SN1992ae* 


- 


- 


-19 


09 


0.06) 





949 


0.027) 


-19 


03 


0.02) 





942 


;o.o2o) 




- 






- 




- 












37.58 


0.390 


E 


(2) 


SN1992al** 


- 


_ 


-19 


42 


0.01) 





984 


U014) 


-19 


28 


0.01) 


1 


030 


;o.oi3) 


-19 


23(0 


02) 


1 


002(0.011) 


-18 


93 




01) 


1 


024 


0.009) 


33.88 


0.292 


- 


(2) 


SN1992aq** 


- 


- 


-18 


95 


'0.03) 





840 


0.034) 


-18 


82 


0.02) 





901 


;0.028) 




_ 








-18 


38 





07) 





818 


'0.050) 


38.25 


0.425 


- 


(2) 


SN1992au* 


- 


- 


-18 


98 


'0.20) 





785 


0.061) 


-18 


91 


0.19) 





844 


'0.134) 




- 






- 


-18 


50 


;o 


03) 





798 


'0.025) 


37.15 


0.487 


E 


(2) 


SN1992bc" 


- 


- 


-19 


60 


'0.03) 


1 


081 


;o.oi5) 


-19 


47 


0.01) 


1 


039 


;0.034) 


-19 


41(0 


03) 


1 


012(0.008) 


-18 


99 


;o 


04) 


1 


027 


0.017) 


34.66 


2.307 


- 


(2) 


SN1992bg** 


- 


- 


-19 


21 


'0.04) 





997 


0.033) 


-19 


10 


0.02) 


1 


025 


;o.oi4) 




_ 








-18 


78 





03) 





981 


0.013) 


35.89 


0.677 




(2) 


SN1992bh 


- 


- 


-18 


84 


'0.02) 


1 


024 


0.022) 


-18 


82 


0.02) 


1 


044 


;0.022) 




- 






_ 


-18 


58 





02) 


1 


028 


0.020) 


36.44 


0.402 


- 


(2) 


SN1992bk* 


- 


- 


-18 


89 


'0.04) 





796 


0.029) 


-18 


83 


0.05) 





806 


;0.025) 




- 






- 


-18 


64 





08) 





740 


'0.045) 


37.00 


0.941 


E 


(2) 


SN1992M* 


- 


- 


-19 


03 


'0.03) 





806 


0.020) 


-18 


94 


0.02) 





836 


'0.013) 




- 






- 


-18 


62 





02) 





806 


0.013) 


36.34 


1.117 


E 


(2) 


SN1992bo* 


- 


- 


-18 


74 


0.02) 





773 


0.011) 


-18 


69 


0.02) 





829 


0.011) 


-18 


75(0 


01) 





781(0.006) 


-18 


52 





01) 





780 


0.005) 


34.49 


0.761 


E 


(2) 


SN1992bp 


- 
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List of parameters of 108 selected SNe la using in this paper. * means "BV bluer" SNe la and ** means "BV bluest" SNe la (see §4.2). SNe la hosted by E or SO galaxies are listed as 
"E". References arc as follows: (1) Wells et al.(1994), (2) Hamuy et al.(1996), (3) Krisciunas et al.(2004b), (4) Altavilla ct al.(2004), (5) Riess et al.(1999), (6) Richmond ct al.(1995), 
(7) Ricss ct al.(2005), (8) Jha ct al.(2006), (9) Krisciunas ct al.(2000), (10) Krisciunas ct al.(2001), (11) Valcntini et al.(2003), (12) Krisciunas ct al.(2004a), (13) Vinko et al.(2003), 
(14) Krisciunas ct al.(2003), (15) VSNET, (16) Kuntal ct al.(2005) 
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SN1999ef** 


-19.04 


; o.io) 


1.423 


0.123) 


-19 


00 





03) 


1 


037 





021) 


-18 


86 





02) 


1 


062 





014) 


-18 


91 





02) 


1 


024 


0.023) 


-18 


49 





03) 





954 





023) 


36.05 


1 


270 




(8) 


SN1999ej* 


-18.68 


0.08) 


0.845 


0.051) 


-18 


33 





03) 





792 





021) 


-18 


26 





01) 





863 





026) 


-18 


32 





02) 





838 


0.014) 


-18 


15 


;o 


02) 





781 


;o 


020) 


33.68 





527 


E 


(8) 


SN1999ek 










-18 


89 





01) 





954 





010) 


-18 


92 





01) 





990 





012) 


-18 


99 





01) 





935 


;o.oi4) 


-18 


82 





01) 





892 





006) 


34.38 





837 




(10) 


SN1999gd 


-17.54 


0.07) 


1.123 


0.069) 


-17 


64 





02) 





942 





010) 


-17 


97 





01) 





984 





010) 


-18 


23 





02) 





956 


0.014) 


-18 


24 





02) 





948 





013) 


34.57 


1 


518 




(8) 


SN1999gh 


-19.01 


0.12) 


0.773 


0.049) 


-18 


62 





02) 





756 





016) 


-19 


12 





03) 





595 





008) 


-19 


07 





01) 





659 


0.005) 


-18 


82 





01) 





692 





004) 


32.86 





946 


E 


(8) 



Table [2] continued. 
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4 RESULTS 



Figure [3] shows the size of typical residuals from the tem- 
plate light curves in U, B, V, R and I-band. Most of ob- 
served light curves follow the templates closely. However, 
note that some early phase U-band photometry (day < 0) 
is systematically brighter than the template. This could be 
due to K-corrections (we made the templates by combining 
observed light curves without K-corrections, see §3.4). 

However, there were a few SNe la which deviated sig- 
nificantly from the templates, due to an intrinsic peculiarity 
in their light curve shapes or a failure to estimate their pho- 
tometric errors properly. Peculiar light curves which don't 
follow the templates are very interesting for studies of the 
nature of supernovae, but that is not our purpose in this 
paper. Therefore, we exclude 14 / 122 SNe la with reduced 
X 2 values larger than 3 (see also Figure [3} from further anal- 
ysis. We did not exclude SNe la due to any other properties 
such as spectroscopic subtypes or host galaxy types. Only 
the exception is that we discuss those peculiar SNe in §4.3 
related to I-band light curves. We refer to the remaining 108 
/ 122 SNe la in the following discussions All parameters of 
this subset of 108 SNe la are given at Tabled As we wrote 
in §2, all of 108 SNe have B and V-band photometry, but 
they don't always have U, R and I-band photometry (see 
Table HJ. 

We also make some subsets of SNe la based on colours, 
stretch factors, host type and redshift in the following dis- 
cussion. Those subsets are all made from the selected 108 
SNe la. 



4.1 Stretch Factor 

In Figure(5] we show distributions of U, B, V, R and I-band 
stretch factors. The values in V, R and I-band are more 
concentrated than those in U and B-band. Since the U-band 
stretch factor suffers from larger errors than the B-band fac- 
tor, due to much fewer measurements and larger photomet- 
ric errors, we conclude that the B-band stretch factor is the 
best choice to serve as an index of SN la properties; we re- 
fer to it extensively in the following sections. Actually, we 
use the reciprocal of the B-band stretch factor in this pa- 
per in order more easily to compare our results with those 
of previous authors who used Amis as a light curve shape 
parameter. 

In Figure [6] we compare our stretch factors with those 
of ALT04, using the common set of 58 SNe la. We obtain 

'WAUvoUla = (°- 86 ± °- 07 ) X V 1 ) + (°- 13 ± °- 07 )- The diS - 
agreement between our stretch factors and those of ALT04 

may be due to the difference in time coverage of templates 
used in the light curve fitting. We find the relationship be- 
tween Amis(B) and the inverse of our B-band stretch fac- 
tor to be Amis = (1-67 ± 0.12) x (sTl - 1) + (1.13 ± 0.02). 
This result is consistent with ALT04, which found Amis = 
(1.98 ± 0.16) x (s^ )AUavilla - 1) + (1.13 ± 0.02). 

Figure [7] shows the relationship between the inverse B- 
band stretch factor and the inverse U, V, R and I-band 
stretch factors. The best fit relations are listed in Table 
[3] JHA06 showed that relations among U, B and V-band 
stretch factors are consistent with a "universal" stretch, and 
our results agree: our B, V, R and I-band stretch factors are 
equal within the formal uncertainties. However, the coefii- 
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Table 4. Average absolute magnitudes (r.m.s.) and stretch factors 
(r.m.s.) 





Magnitude 


Stretch Factor 


Number 


U-band 


-19.08 (0.79) 


1.07 (0.14) 


44 


B-band 


-18.79 (0.76) 


0.96 (0.12) 


108 


V-band 


-18.84 (0.58) 


1.00 (0.13) 


108 


R-band 


-18.93 (0.53) 


1.00 (0.13) 


88 


I-band 


-18.70 (0.45) 


0.98 (0.15) 


102 



cient of the relationship between U and B-band stretch fac- 
tors is smaller than those for the relationships between B, V, 
R and I-band stretch factors. If we select 19 SNe la in com- 
mon with JHA06, the relation between U and B-band stretch 
factors is s (u) = (1.03+0.11) x s (B) +(0.05+0.11). This result 
is consistent with those of JHA06 (s(in = 1.04 x S( fl ) —0.02). 



4.2 Magnitude and colour 

Figure [S] shows the distribution of U, B, V, R and I-band 
absolute pea k magnitudes c o rrecte d for Galactic extinction 
according to ISchleeel et all ljl998h . The trend that redder 
bands are more concentrated than bluer bands possibly in- 
dicates the effects of extinction by dust (see also Table [4}. 

The colours of the SNe la also suggest that dust plays 
an important role. Figure shows the distribution of the 
(B — V) colour at the time of B-band maximum, (B — V) ma x- 
The tail of redder side in Figure[9]are likely SNe la reddened 
by dust in their host galaxies. As shown in Figure [TD] redder 
SNe la tend to be missing at higher redshift, implying that 
the redder SNe la are affected by dust. The clear edges of the 
blue side of the distributions in Figures and [TU] indicates 
that the (B — V) m ax colour of bluest SNe la is almost con- 
stant. Fitting a single-sided gaussian to the left side of the 
histogram in Figure — V) m ax > —0.20), we estimate 
the intrinsic colour of bluest SN la as (B — V) m ax = —0.12. 
Figure [TT] shows relationships between (B — V)max and the 
other colours at the time of B-band maximum. We will dis- 
cuss this in detail in §5. 

We should pay attention to the possibility that there 
are intrinsically redder SNe la. We cannot distinguish those 
redder SNe la from the SNe la reddened by host galaxy dust 
based on apparent colour. In §6 and §7, we discuss about the 
issue. 

Based on the the result of single-side gaussian fitting of 
B-V distribution, we made three subsamples, "BV bluest" 
SNe la (-0.14 < (B — V) max «S -0.10, -0.12 ± 0.02 
mag around the median), "BV bluer" SNe la (—0.10 < 
(B — V)max < —0.02, 0.1 mag redder than the median) 
as less dust affected samples and "BV redder" SNe la 
(0.00 < (B — V) m ax) as a dust affected sample. We use 
them in the following discussions. 



4.3 I-band Peculiarity 

It is well known that there is quite a variety of shapes in the 
I-band light curves of SNe la. Usually, there is a "bump" or 
a "second peak" in the I-band, but some SNe la don't have 
a secondary peak (see Figure [12]). Some of previous authors 



fitted these I-band light curves using two Gaussian functions 
jContardo et al.ll2000l : iNobili et al.ll2005l ). However, we tried 
to fit I-band light curves using only one template to avoid 
introducing extra parameters. 

From all of 122 SNe la, we selected 71 well observed SNe 
la which have at least one I-band measurement in two peri- 
ods, -10 days to days and +10 days to +20 days. Based on 
photometry during these periods, the 71 SNe la were classi- 
fied into two types, normal I-band light curves and peculiar 
I-band light curves (see Figure [T3")) . where "peculiar" means 
that the I-band light curves are fainter than the I-band tem- 
plate from -10 days to days and brighter than the I-band 
template from +10 days to +20 days. 

We found about 90% of SNe la (63/71 SNe) can 
be fitted well with the template and only 8 SNe la 
(SN1991bg, SN1995ac, SN1997br, SN1998ab, SN1999ac, 
SN1999by, SN1999gh, SN2002cx) can not be fitted well with 
the template. Most of those which deviate from the template 
are known to be spectroscopically peculiar as well. For exam- 
ple, SN1991bg and SN1999by are classified as "SN1991bg- 
like," SN1995ac and SN1997br are classified as "SN1991T- 
like," and SN 2002cx is known as very peculiar SNe la 
|Li et al.ll2003l ). On the other hand, I-band light curves of 
SN1991T and SN1995bd, which are classified spectroscop- 
ically as "SN1991T-like" do fit the I-band light template. 
Based on the value of reduced \ 2 calculated with Multi-band 
Stretch method, these SNe which do not fit the template I- 
band light curve have been excluded from the samples in 
the follo wing discussions, as we mentioned in the beginning 
of §4. As lNobili etail (|2005l ) pointed out the correlation be- 
tween lack of a second I-band maximum and a lower B-band 
stretch, the Multi-band Stretch Method would have room for 
improvement in the treatment of these SNe. 



5 SUB GROUP OF TYPE IA SUPERNOVAE 

Do all SNe la with the same B-band stretch factor have the 
same intrinsic colours? The best way to answer this question 
is to compare a series of well-calibrated spectra from differ- 
ent events, but we lack such a dataset. Therefore, we have 
investigated the issue using photometric information from 
our sample of SNe la. 

In Figure 1141 we plot SNe la classified on the basis of 
their (B — V) ma x colour as "BV bluest" SNe la, "BV redder" 
SNe la and "others". As we discussed in §4, the bluer SNe 
la may be free from dust extinction. In Figure [141 there is a 
"sub" group at the lower right, with inverse B-band stretch 
factor> 1.1 and (V— I) ma x < 0. This "sub" group consists of 
SNe la which are redder in (B—V) ma x, though they are bluer 
in (V — I)max- It is difficult to explain this combination of 
colours using ordinary SNe la shining through dust in their 
host galaxies. Instead, it may mean that the "sub" group is 
due to a set of SNe la with intrinsic SEDs which differ from 
the "main" group with inverse B-band stretch factor< 1.1, 
at least in B, V and I-band. 

There are some observational indications that proper- 
ties of SNe la are related to their environment. lHamuv et all 
jl996t ) showed that the mean luminosity of SNe la in ellip- 
tica l galaxies is fain t er tha n tha t of SNe la in s piral g alax- 
ies. iMannucci et all (|2005h and ISullivan et all (|2006l ) also 
showed that the star-forming activity in a galaxy is related 
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Table 3. Relationships between inverse B-band stretch factor and 
stretch factors in U,V, R and I-band 



Band Relation r.m.s. (mag) Number 



U-band 


(0.78 ±0.07) 


X 


S (_B) 


+ (0.14 ±0.07) 


0.040 


44 


V-band 


(0.93 ±0.05) 


X 


S (B) 


± (0.03 ± 0.05) 


0.009 


108 


R-band 


(1.04 ±0.05) 


X 


S (B) 


- (0.05 ± 0.05) 


0.009 


88 


I-band 


(1.02 ±0.05) 


X 


S (B) 


- (0.01 ±0.06) 


0.011 


102 



to properties of its SNe la. We have plotted the SNe la 
hosted by E or SO galaxies in Figure 1151 We note that the 
SNe la of the "sub" group at the lower right tend to be 
found in E or SO galaxies. However, not all SNe la found in 
E or SO galaxies belong to this group. 

Of course, a difference of intrinsic colours doesn't mean 
a difference of intrinsic luminosity directly. However, a dif- 
ference in intrinsic colours affects the accuracy of correcting 
for host galaxy dust extinction if we correct the extinction 
according to apparent colour of SN la. For example, the 
bluer SNe la with inverse B-band stretch factors > 1.2 in 
the (B — V) panel of Figure [151 are redder than the bluer SN 
la with inverse B-band stretch factors < 1.2 by about 0.05 
magnitude in (B — V) ma x- If we treat the (B — V) m ax as a 
constant when we estimate the host galaxy dust extinction, 
the correction will be overestimated systematically by about 
0.2 magnitude in B-band. This is critical if we want to use 
SNe la as distance indicators for cosmological studies. 



6 RELATIONS BETWEEN B-BAND STRETCH 
FACTOR AND BRIGHTNESS 

Figure [TH] and Figure [T7] show the relationships between in- 
verse B-band Stretch factor and luminosity at the B-band 
maximum (hereafter, stretch - magnitude relation) in U, B, 
V, R and I-band. The scatter in the luminosity is mainly 
caused by three factors: the peculiar velocity of the host 
galaxy, extinction by dust in the host galaxy, and intrinsic 
diversity of SNe la. 

Peculiar velocity makes the dispersion larger when we 
estimate the distance to a SN from its host's recession veloc- 
ity. The effect of peculiar velocity must be treated properly 
if the host galaxy is nearby. For example, if a host galaxy has 
peculiar velocity of 400 km/s, we will overestimate or under- 
estimate the luminosity of its SN by ~0.3 mag at z = 0.01 
and ~0.1 mag at z = 0.02. However, statistically speaking, 
the effect can be small if the sample size is large enough. 
In our case, Table [S] shows that there is no significant ef- 
fect on averaged absolute magnitudes caused by peculiar 
velocity if we choose "BV bluest" + "BV bluer" SNe la in 
order to avoid the effect of host galaxy dust extinction. In 
Figure 1181 we plotted individual SNe la shown in Table [5j 
however there are no considerable effect caused by peculiar 
velocities. Table \S\ also shows the scatter size in the nearest 
sample (0.00 < z ^ 0.02) is only 10% larger than the further 
samples. 

It is better to include properly the effect of peculiar 
velocity in error bars, however we don't know which value 
could represent peculiar velocities. For example, peculiar ve- 



Table 5. Average magnitudes of "BV bluest" + "BV bluer" SNe la 
at different redshifts 



Redshift 


M B 


Number 


0.00 < z < 0.02 


-19.13 (0.35) 


18 


0.02 < z < 0.05 


-19.15 (0.23) 


19 


0.05 < z 


-19.01 (0.26) 


10 



locity of galaxies in the clusters often reaches ~l,000km/s. 
It is very difficult to estimate appropriate value of pecu- 
liar velocities without additional observations, and we don't 
want to unnecessarily overestimate the size of error bars. 

One way to avoid the peculiar velocity problem is to 
use further redshift SNe la only. But, the number of SNe la 
becomes smaller in compensation for the selection. Hence, 
we didn't make any corrections for the effect of peculiar 
velocity in this study. 

Extinction due to dust in the host galaxy is a major 
contributor to the scatter in luminosity. Unfortunately, it is 
very difficult to estimate this extinction properly. We must 
separate the effect of host galaxy dust extinction from in- 
trinsic colour variations of SNe la. One way to avoid this 
problem is to use only bluer SNe la which are expected 
to be free from extinction. In Figure 1161 we plot the "BV 
bluest" SNe la and "BV bluer" SNe la samples. Coefficients 
of the relationship between absolute magnitude and inverse 
stretch factor are listed in Table [6] In all passbands, the 
slope of the line derived from the "BV bluest" sample is 
about double that of the slope of the line derived from the 
combined "BV bluest" + "BV bluer" sample. Figure [TBI and 
Figure [17] also show that most of "BV bluest" SNe have 
broader light curve widths than "BV bluer" SNe. This sug- 
gests there may be different groups of SNe la with different 
intrinsic properties. The slope of the stretch - magnitude re- 
lation of ALT04 is 1.83 (using 1.102 in Table 1 of Altavilla's 
paper and Amis = (1.67±0.12) x (s7^ - 1) + (1.13±0.02)), 
which is shallower than the slope of "BV bluest" SNe la but 
steeper than the slope of the combined "BV bluest" + "BV 
bluer" sample. 

We also tried to compare the dispersion of luminosity of 
"BV bluest" SNe la around the stretch - magnitude relation 
derived from "BV bluest" SNe la with that of SNe la hosted 
by E or SO galaxies. Those SNe la hosted by E or SO galaxies 
should be largely free from extinction due to interstellar dust 
in the host galaxy, but Figure [T5] shows that they have a 
larger dispersion than "BV bluest" SNe la (which are also 
expected to be free from host galaxy dust extinction); see 
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also Table [7] Even if we select "BV bluer" SNe la hosted 
by E or SO galaxies, the dispersion is still larger than "BV 
bluest" SNe la without any dust / colour corrections. These 
facts may indicate that SNe la hosted by E or SO galaxies 
have larger intrinsic diversity than SNe la hosted by other 
type galaxies, or that the host galaxies have some significant 
dust extinction. 

Of course, the scatter of brightness of SNe la in E or 
SO galaxies could be attributed to the peculiar velocity of 
the host galaxies, but the peculiar velocity of host galaxies 
can not account for the fact that direction of distribution 
of SNe la hosted by E or SO galaxies shown in Figure [TT] is 
parallel to the direction of distribution due to dust in the 
Milky Way. 

6.1 Colour Correction For Absolute Magnitude 

It is known empirically that correcting the luminosity of SNe 
la based on their colours will decrease the dispersion of peak 
B-band luminosities. For example, MLCS, CMAGIC, and 
SALT use colour information as part of the light curve fitting 
technique, while Stretch and Amis methods use colour infor- 
mation to decrease the dispersion after fitting light curves. 
We also introduce colour terms to decrease the dispersion 
around the stretch - magnitude relation. 

We need a colour excess to remove the extinction due to 
dust in the host galaxy, and the intrinsic colour of each SN 
la to correct for diversity among SNe la. However, it is diffi- 
cult to separate an observed colour excess into the intrinsic 
colour of the SN and the reddening by dust. Therefore, we 
don't try to distinguish them; instead, we use the observed 
colour of SNe la empirically to reduce the dispersion of the 
relationship between absolute magnitude and inverse stretch 
factor. We employ a chi-square fitting technique with follow- 
ing equation. 

X = M B - {ax s^ 1 , + 6 x {B — V) max + c} (3) 

Here Mb is the absolute B-band magnitude at maxi- 
mum shown in Table [2] without any corrections, (B — V) ma x 
is the apparent colour at B-band maximum and a, b, c are 
fitting parameters. Note that here we treat (B — V) m ax as 
independent from B-band stretch factor. 

Table [8] shows the inverse B-band stretch factor and B- 
band luminosity relationships, including the r.m.s. of each 
sample. The r.m.s. is the smallest in sample 8D, which in- 
cludes SNe la hosted by E or SO galaxies at higher redshift 
(z > 0.02). The relationship of sample 8D is similar to that 
of sample 8B (SNe at z > 0.02), but the r.m.s. is smaller in 
the sample drawn from E or SO galaxies. This may indicate 
that extinction by dust in an E or SO host galaxy is smaller 
than extinction in other types of galaxies, if both sample 8B 
and 8D have the same intrinsic dispersion of peak B-band 
luminosity. 



7 HOST GALAXY DUST PROPERTIES 

There are several approaches to estimate the effects of ex- 
tinction by dust in host galaxies. In this section, we esti- 
mate host galaxy dust extinction using the apparent colours 
of SNe la. Dust in a host galaxy dust makes the apparent 
colour of a SN la redder, so we can estimate how a host 



Table 7. r.m.s. of the stretch - magnitude relation of the samples and 
their average (B — V) ■max 



Sample 


r.m.s. (mag) 


average (B - V) m ax 


Number 


7A 


0.17 


-0.12 


21 


7B 


0.51 


0.04 


27 


7C 


0.30 


-0.07 


15 


7D 


0.17 


-0.12 


3 



7A is "BV bluest" sample. 

7B is sample hosted by E or SO galaxies with no colour selection. 
7C is "BV bluer" sample hosted by E or SO galaxies. 
7D is "BV bluest" sample hosted by E or SO galaxies. 



galaxy obscures a SN using colour information. We use the 
following relationship: 

A x = RxX E(M B - M v ) (4) 

Here Ax is dust absorption in X-band, Rx is the ratio 
of total to selective extinction in X-band, and E(Mb — My) 
is the colour excess. We derive E(Mb — Mv) assuming an 
intrinsic (B — V) ma x colour for each SN la and extract the 
excess from observed colour. 

There are two main problems when we use equation (4) 
to correct for extinction in the host galaxy. One problem is 
that we don't know if the properties of dust in host galaxies 
are the same as those of dust in the Milky Way, especially 
the value of Rx- ALT04 and other works showed that Rb 
of dust in host galaxies may be smaller than Rb of dust in 
the Milky Way. Another problem is that we don't know the 
intrinsic colour of SNe la. Usually we assume that the bluest 
SNe la are free from extinction, but there is no guarantee 
that all SNe la have the same intrinsic colour as the bluest 
examples. These two problems make the correction difficult. 
Hence, we reduce the problem to two extreme cases in the 
following. 

In case A, we assume that the intrinsic colour of SNe 
la depends on the inverse B-band stretch factor. This as- 
sumption was also used in ALT04. We found earlier that 
the intrinsic (B — V) m ax = —0.30 x s™ + 0.18 from the 
relations between and B, V-band (see Table [BJ. This 
assumption is consistent with SNe la with s7^ < 1.1, but 
not consistent with SNe la with > 1.1 (see Figure [15)) . 

In case B, we assume the intrinsic colour of SNe la is 
constant, and has some significant scatter. In this case, even 
if two SNe la have the same inverse B-band stretch factor, 
they don't necessarily have the same colour. 

First, we tried to correct for extinction due to dust in 
host galaxies under the assumption of case A, using three dif- 
ferent values for Rb '- 4.315 (appropriate for the Milky Way, 
ISchlegel et allll998l l. 3.5 (ALT04) and 3.0. We selected red- 
der SNe la, which are (B — V) m ax > —0.02 (0.1 mag redder 
from average "intrinsic" (B — V) m ax = —0.12 derived from 
Figure [9} • These SNe la must be significantly reddened by 
dust in the host galaxy. Table [9] shows the scatter around 
the stretch - magnitude relationship in B-band luminosity 
after we have corrected for extinction according to equation 
(4). Though smaller values of Rb do decrease the disper- 
sion, the scatter is still large. This means we must introduce 
another factor to account for the large dispersion; for ex- 
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Table 6. Relationships between B-band stretch factor and luminosity 



Sample Relation r.m.s. (mag) Number 



U-band 


6A 


(3.28 ±0.40) x 


'W 


- (22.93 ± 0.42) 


0.32 


8 




6B 


(1.50 ±0.40) x 




- (21.11 ±0.44) 


0.32 


18 


B-band 


6A 


(2.28 ±0.42) x 


4) 


- (21.49 ±0.41) 


0.17 


21 




6B 


(1.08 ±0.24) x 


s w 


- (20.30 ± 0.24) 


0.23 


46 


V-band 


6A 


(2.58 ±0.38) x 




- (21.67 ±0.38) 


0.20 


21 




6B 


(1.10 ±0.24) x 




- (20.22 ±0.25) 


0.22 


46 


R-band 


6A 


(1.88 ±0.57) x 




- (20.96 ±0.56) 


0.18 


14 




6B 


(0.75 ±0.29) x 




- (19.89 ±0.30) 


0.26 


34 


I- band 


6A 


(1.10 ±0.61) x 


s w 


- (19.93 ±0.60) 


0.23 


20 




6B 


(0.52 ±0.30) x 


v 1 ) 


- (19.41 ±0.31) 


0.27 


42 



6A is "BV bluest" sample. 

6B is "BV bluest" +"BV bluer" sample. 



Table 8. Relationships between B-band inverse stretch factor and luminosity with colour 
term and r.m.s. of the stretch - magnitude relation 



Sample relation r.m.s. (mag) Number 



8A 


(0.96 ± 0.28) 


x 


S (_B) 


- (2.51 ±0.15) 


X 


(B- 


V)max 


- (20.26 ± 0.29) 


0.48 


108 


8B 


(0.98 ±0.18) 


x 


S (B) 


- (2.28 ±0.28) 


X 


(B- 




- (19.95 ± 0.19) 


0.27 


50 


8C 


(1.09 ±0.59) 


x 


S (S) 


- (1.78 ±0.43) 


X 


(B- 




- (20.15 ± 0.71) 


0.33 


27 


8D 


(0.99 ±0.34) 


x 


S (B) 


- (2.23 ±0.24) 


X 


(B — 




- (20.10 ± 0.41) 


0.12 


15 



8A is all SNe la. 

8B is SNe la of z > 0.02. 

8C is SNe la hosted by E or SO galaxies. 

8D is SNe la of z > 0.02 hosted by E or SO galaxies. 



ample, some intrinsic scatter of (B — V) max . Note that the 
r.m.s. of all samples after our corrections for dust in the host 
galaxy is larger than the r.m.s. shown in ALT04, Figure 5 
(b) (Rb = 3.5, a = 0.31). The explanation is mainly that 
we used only redder SNe in our analysis, but the ALT04 
sample includes bluer SNe la which have smaller scatter in 
B-band luminosity than redder SN la. If we use a set of 55 
SNe la in common with ALT04 to derive the dispersion, we 
find a = 0.44 mag with Rb = 0.35, which is smaller than 
that of redder SN la (R B = 3.5, a = 0.81, Table but 
still larger than the result shown in ALT04, Figure 5 (b) 
(Rb = 3.5, a = 0.31). The reason may be in the details of 
how one derives E(B-V): ALT04 used a weighted average of 
(B - V)max and (B - V)tail, but we used (B - V) max only. 

Next, we tried to correct host galaxy dust extinction 
under the assumption of case B (a constant intrinsic colour 
for all SNe la) using a factor Rb = 4.315. Indeed, according 
to Figure [TT1 there is no basis to assume that dust in host 
galaxies has properties different from dust in the Milky Way. 
Under the assumption of case B, we can't derive Ax from 
equation (4) because we don't know the intrinsic (B — V) m ax 
of SNe la. Instead, we tried to derive E(Mb — My) from the 
following equation (5) , under the assumption that the intrin- 
sic B-band luminosity of SNe la is given by the relationships 
listed in Tabled 



Table 9. Dispersion of B-band luminosities (mag) of selected 
redder SNe la after host galaxy dust correction based on three 
different conversion factors of Rb ■ 



Sample 


R B = 4.315 


R B = 3.5 


Rb = 3.0 


Number 


all 


0.89 


0.81 


0.79 


53 


9A 


0.99 


0.90 


0.89 


38 


9B 


0.56 


0.48 


0.45 


15 


ALT04 


0.52 


0.44 


0.42 


55 



9A is "broader" LC sample, which s^b) 

9B is "narrower" LC sample, which s7^. > 1.1. 

ALT04 is common 55 samples with ALT04, Figure 5 (b). 



M B - (2.28 X s7± - 21.49) 

E(M B - My) = — -J2 L ( 5 ) 

Hmw 

Here we regard residuals from the intrinsic B-band lu- 
minosity at some particular B-band stretch factor as due 
completely to extinction by dust in the host galaxy, and we 
divide the residuals by Rmw = 4.315, the ratio of total to 
selective extinction for dust in the Milky Way. We applied 
equation (5) to dust-affected SNe la, corrected them for red- 
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dening, then compared their (B — V) ma x distribution with 
that of dust-free SNe la. 

To make the dust-affected sample, we selected SNe la 
which are fainter by at least 1.0 magnitude in B-band than 
the intrinsic luminosity derived from inverse B-band stretch 
factor. For the dust-free sample, we selected SNe la which 
fall within 0.17 magnitude (which is one standard deviation 
of "BV bluest" SNe la) of the relationship between intrinsic 
luminosity and inverse B-band stretch factor. 

The results are distributions of (B — V) colours for 

V / max 

dust-free samples, dust-affected samples, and dust-corrected 
samples, using B-band residuals to correct for extinction 
(Figure I20p . or V-band residuals to correct for extinction 
(Figure [2lT) . We compare the distribution of SNe la which 
are corrected using the factor Rb of the Milky Way to the 
distribution of SNe la which are corrected with a smaller 
factor Rb- These comparisons show that if we accept the 
idea that intrinsic SNe la colours have an intrinsic dispersion 
which is similar to that of the "reference" samples in Figure 
1201 and 1211 the average factor R of dust in host galaxies is 
not different than that for dust in the Milky Way. We ran 
a Kolmogorov-Smirnov test on the distribution of corrected 
(B — V) colours and found no significant difference between 
them. Intrinsic diversity among SNe la diver sity has also 
been s hown by previous spectroscopic studies. iBerian et al.l 
(2006), for example, showed that r.m.s. of B-band and V- 
band luminosities is about 0.1 - 0.2 magnitude, and their 
results are consistent with ours. 

If the intrinsic colours of SNe la have a significantly 
large scatter, it is very difficult to correct properly for ex- 
tinction by dust in host galaxies with photometric measure- 
ments only. 



magnitude relationship of colour-corrected SNe la hosted by 
distant E or SO galaxies (8D, Table[8]) is 0.12 mag, which is 
smaller than that of the "BV bluest" SNe la discussed above. 
However, the properties of SNe la in E or SO galaxies tend 
to differ from those of SNe la within other types of galaxies 
(cf. Figure [T?] Figure fUJl and Figure [T5)l . This means we 
may have to apply different relations to SNe la hosted by E 
or SO galaxies. 

Even if we select SNe la within early type galaxies, those 
SNe are not always free from dust, as shown in Figure 1111 
This is consistent with the idea that not only interstellar 
dust but al so dust s urrou nding the SN a ffects the light from 
SNe la, as IWang I (|2005l 1 pointed out. iPatat et al'l (|2007T l 
reported the spectroscopic detection of circumstellar dust 
around normal SN la. The SNe la occuring in early type 
galaxies also tend to have narrower light curve widths, with 
inverse B-band stretch factor > 1.1 (see Figure [15}, and 
they have larger intrinsic dispersion in SED as shown in Fig- 
ure [15] Figure [191 shows these SNe la in early- type galaxies 
have larger dispersion than "BV bluest" SNe la, or that the 
stretch - magnitude relation is different from that of "BV 
bluest" SNe la. Our results are consistent with the distri- 
bution of B-band stretch f actors of nearby SNe la within 
E or SO hosts determ ined in lSullivan et al.l ()2006h . However, 
ISullivan et all (|2003l ) also claimed that at high redshifts, the 
scatter in luminosity of SNe la hosted by E or SO galaxies is 
smaller than that of SNe la inside spiral galaxies. 

In any case, we must know whether or not the proper- 
ties of SNe la in the low-z universe and the high-z universe 
are the same, since it is important when we apply stretch- 
brightness relationships or make K-corrections. To answer 
these questions, we need to study the nature of SNe la at 
all redshifts. 



8 SN IA SAMPLES FOR COSMOLOGY 

In §6, we studied the connection between light curve shape 
and luminosity and looked for ways to decrease the disper- 
sion around the relationship. In §7, we checked the properties 
of dust in host galaxies. 

These results suggest that if one wishes to use SN la 
as standard candles, one should take one of the following 
approaches. First, one might select only "BV bluest" SNe 
la. As we showed in §6.1, colour corrections make the dis- 
persion smaller, but the improvement in the dispersion of 
a combined sample of "BV bluer" and "BV bluest" SNe la 
after colour corrections is only 0.02 mag, and their r.m.s. 
is still larger than that of the "BV bluest" sample with no 
colour correction. In addition, Figure [15] shows that the in- 
trinsic colour of SNe la seems to have a large scatter, and 
that even SNe which have same B-band stretch factor may 
not have same intrinsic colour. Until we solve two puzzles, (i) 
how to separate the apparent colour of SNe la into intrinsic 
colour and reddening due to the host galaxy, (ii) the value of 
the ratio of total to selective extinction R and its evolution 
from the low-z to the high-z universe, we had better not 
use colour corrections in samples which have large colour 
diversity. The study of K-corrections is very important for 
selecting "BV bluest" SN la and determining extinction by 
dust in host galaxies. 

Another approach is to select SNe la which occur in 
E or SO galaxies. As shown in §6.1, the r.m.s. of stretch - 



9 SUMMARY 

In this work, we re-analyze U, B, V, R and I-band light 
curves of 122 previously-published near by (redshift < 0.11) 
SNe la, which are main ly from JHA06, iRiess et al l (|l999aT ) 
and lHamuv et all l|l996Tl . We devise the "Multi-band Stretch 
method", which is a revised Stretch method extended to 
fitting up to five bands. We create new U, B, V, R and I- 
band light curve templates from well-observed SNe la for 
our "Multi-band Stretch method". It is known that the I- 
band light curves of SNe la shows much variation. However, 
we found our I-band template fitted about 90% of SNe la 
I-band light curves very well. 

Using the results of our light curve fitting, we exam- 
ine the relationships between luminosity, colours at maxi- 
mum light, and stretch factors. In §5, we show a possible 
sub-group of SNe la which have different colours from the 
main group of SN la; this sub-group tends to be found in 
E or SO galaxies. In §6, we find the relationship between 
inverse B-band stretch factor and U, B, V, R and I-band 
luminosity, using colour-selected SNe la to avoid contam- 
ination from the outlying sub-group. We obtain a stretch 
- magnitude relationship for "BV bluest" SNe la without 
colour correction as (2.28 ± 0.42) x sjL ~ (21.49 ± 0.41), 
which is steeper than the slope of the corresponding stretch- 
magnitude relation in ALT04. The r.m.s. of the residuals 
from this relationship is 0.17 mag. If we attempt to re- 
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duce the dispersion in luminosity by using colour correc- 
tions, we find that for distant (redshift > 0.02) SNe la, 
Mb = 0.98 x s-g } - 2.28 x (B - V) max - 19.95 with an r.m.s. 
of 0.27 mag. The r.m.s. becomes 0.12 mag when we select 
only SNe la hosted by E or SO galaxies. In §7, we statisti- 
cally investigate the properties of dust in host galaxes, and 
find that the ratio of total to selective extinction R can be 
consistent with that of dust in the Milky Way if we assume 
that SNe la have some diversity in their intrinsic colours. 

Based on these results, in §8, we discuss how to select 
subsets of SNe la which may be good distance indicators for 
cosmology. We find two possibilities: one is to use only "BV 
bluest" SNe la with broad light curves, and the other is to 
use only SNe la inside E or SO galaxies. To use these samples 
for precision cosmology, we need to study the properties of 
SNe la from the low-z to the high-z universe. 
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Figure 1. Redshift distribution of SNe la used in this work. Most 
of them arc at z < 0.02. 
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Figure 3. Residuals from light curve template fitting in U, 
B, V, R and I-band magnitude. All epochs of 122 SNe la are 
plotted after fitting to the templates. Variance index is defined 

aS Am X I \J template + photometry' h ° TC > Am X is thc residual 

between the X-band template and the observation (plus sign 
means that photometry is fainter than template). a- temp i ate and 
^photometry are typically in the order of 0.01 mag. Early phase 
U-band photometry (day < 0) are systematically brighter than 
the template. A possible reason of the excess is K-corrcction. 



Figure 2. U, B, V, R and I-band light curve templates (stretch 
factor = 1) with scatter. The r.m.s. (1<t) are derived from residuals 
of composited light curves and shown in dashed line. The B, V, R, 
I-band templates are from -10 days to +80 days, and the U-band 
template is from -10 days to 70 days after B-band maximum. The 
U-band template has larger dispersion than B, V, R and I-band 
templates because U-band photometry is poorer than B, V R and 
I-band photometries. 
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Figure 4. Distribution of reduced \ 2 °f the fitting. Most of 
SNe la well fit in the templates, but there are some SNe la 
not well fitted in the templates. x 2 / d -°-f- > 4 : SN1994T, 
SN1996bo, SN1997cn, SN1998aq, 



SN1999by, SN1999da, SN2002cx, SN2003du. x 2 /d.o.f. 



SN1998de, SN1999ao, 
> 3 

: SN1992A, SN1992ag, SN1998bp, SN1999aw. A few of them 
are known as spectroscopically peculiar SNe la, but some of the 
others are known as "Branch normal" SNe la. 
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Figure 6. Comparison with ALT04. The left panel shows a re- 
lation between inverse B-band stretch factors. The right panel 
shows the relation between inverse B-band stretch factor (this 
work) and Amis from ALT04. The solid line is a regression line 
derived from least x 2 fitting. 



> 



o 

c 
o 
— 

CD 

u 




0.5 



0.5 1 1.5 0.5 1 1.5 

Stretch Factor(B) -1 



Figure 7. Relations between the inverse B-band stretch factor 
and inverse U, V, R, I-band stretch factors of all SNe la. The 
solid lines are regression lines (see also Tabl J31 . 



Figure 5. Histograms of the stretch factors in U, B, V, R and 
I-band. There are 44 SNe la in U-band, 108 SNe la in B-band, 
108 SNe la in V-band, 88 SNe la in R-band, 102 SNe la in I-band. 
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Figure 8. Histograms of absolute magnitudes in U, B, V, R and I- 
band. The magnitudes arc corrected for Galactic dust extinction, 
but not for host galaxy dust extinction. Bluer bands such as U 
and B-band have broader distributions than redder bands. 




Figure 10. (B — V) m ax distribution as a function of redshift. 
Error bars include only photometric errors and fitting errors (pe- 
culiar velocity of galaxies are not included). Redder SNe la tend 
to be missing at larger redshift because they are fainter than bluer 
SNe la which expected to be free from host galaxy dust. Also this 
figure shows there are no apparent trend in the bluest end of 
(B — V)max colour along redshift at the range of z < 0.1. 



Figure 9. Rest frame distribution of [B — V)max without host 
galaxy dust correction. Tail in redder side can be attributed host 
galaxy dust reddening. We can regard bluer SNe la may be less 
affected by host galaxy dust. Median (B — V) m ax of one-side 
Gaussian fit using only bluest SNe la is -0.12 +/- 0.02. This value 
is consistent with ALT04. 
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Figure 11. Colour-colour diagram of SNe la at B-band maxi- 
mum. Filled red squares are SNe la hosted by E or SO galaxies. 
Arrows show direction and size of reddening when we assume that 
absorption properties of host galaxies dust are same as Galactic 
dust (Ru = 5.434, R B = 4.31 5, R v = 3.315, R R = 2.673, Rj = 
1.940 from lSchlegel et al,lll998l) and Ab = l.Omag. Even if we se- 
lected only SNe la hosted by Elliptical galaxies, they are plotted 
along the Galactic dust reddening direction. The colour have been 
corrected for Galactic dust reddening, but not for host galaxy dust 
reddening. 



Figure 12. Examples of normal I-band light curve and peculiar 
I-band light curve. SN1994ae is a spectroscopically normal SN la, 
and there is a "bump" in the I-band light curve. SN1999by is a 
spectroscopically peculiar SN la, and there is no "bump" in the 
I-band light curve. 
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Figure 13. I-band light curves of 71 selected SNe la (top panel) 
and distribution of the epochs (bottom panel). Variance index 

is defined as Am//y/^^ + T 2 photometry , where Am, is 
the residual between the I-band template and the observation, 
^photometry is the photometric error and cr ternvlate is the tem- 
plate error. We classified the SNe into two types, normal I-band 
light curves (open square) and peculiar I-band light curves (filled 
red circle), according to photometry points (see text). 
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Figure 14. Distribution of (V — I) m ax as a function of the inverse 
B-band stretch factor. Filled blue circles are "BV bluest" SNe la, 
red crosses are "BV redder" SNe la and open triangles are the 
other SNe la. 



Figure 15. Distribution of rest frame colours as a function of the 
inverse B-band stretch factor. Filled red squares are SNe hosted 
by E or SO galaxies. The colour have been corrected for Galactic 
dust reddening, but not for host galaxy dust reddening. 



This paper has been typeset from a TjrjX/ J^TgX file prepared 
by the author. 
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Figure 16. Relations between inverse B-band stretch factor and 
absolute B-band magnitude. Filled blue circles are "BV bluest" 
SNe la, open blue squares are "BV bluer" SNe la and the other 
symbols are for other SNe la. The solid line is a relation between 
inverse B-band stretch factor and B-band magnitude derived from 
"BV bluest" SNe la, while the dashed line is a relation derived 
from "BV bluest" + "BV bluer" SNe la. Redder and fainter SNe 
la than those " B V bluest" or " B V bluer" SNe la may be affected 
by host galaxy dust. 



Figure 17. Relations between inverse B-band stretch factor and 
absolute magnitude of U, V, R and I-band. Filled blue circles 
are "BV bluest" SN la, open blue squares are "BV bluer" SNe 
la and the other symbols are for other SNe la. The solid line 
is a relation between inverse B-band stretch factor and B-band 
magnitude derived from "BV bluest" SNe la, while the dashed 
line is a relation derived from "BV bluest" + "BV bluer" SNe la. 
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Figure 18. Absolute magnitude distribution as a function of red- 
shift. Magnitude is corrected for Galactic dust extinction, but not 
for host galaxy dust extinction. The bright outlier near z = 0.007 
is SN1991T. If peculiar velocity is not negligible, the number of 
brighter SNe la than the average increase in lower-z (also fainter 
SNe la increase, but we can't distinguish them from the SNe la 
which extincted by host galaxy dust). However, there is no con- 
siderable effect caused by peculiar velocity. See also Table [5] 



Figure 19. Residuals from the stretch - magnitude relation de- 
rived from "BV bluest" SNe la without dust /colour correction. 
Filled blue circles are "BV bluest" SNe la (sample 7A in Table 
0, red symbols are SNe la hosted by E or SO galaxies (sample 
7B in Table 0, filled red triangl es are "BV bluer" SNe la hosted 
by E or SO galaxies (sample 7C in Table 0, open red triangles 
with filled blue circle are "BV bluest" SNe la hosted by E or SO 
galaxies (sample 7D in Table and the other symbols are for 
other SNe la. 
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Figure 20. Distributions of (B-V) m ax of dust free and dust cor- 
rected SNe la using B-band brightness. The top panel is reference 
SNe la which are mostly dust free. The second panel is for SNe la 
which are significantly affected by host galaxy dust, but without 
dust correcti ons. The third pane l is after dust correction using 
Rb = 4.315 llSchlegel et alj|l99Sl) . and the bottom panel is dust 
correction with Rb = 3.5 (ALT04). Shaded bars show type la 
SNe with broader light curve shape, < 1.1, which are ex- 

pected to be spectroscopically "normal" SNe la (see Figure [14J. 



Figure 21. Distribution of (B-V)max of dust free and dust cor- 
rected SNe la using V-band brightness. The top panel is reference 
SNe la which are mostly dust free. The middle panel is for SNe la 
which are significantly affected by dust, but without dust correc- 
tions. The bottom p anel is after dust correction with Ry = 3.315 
llSchlegel et al.ll998h . Shaded bars show SNe la with broader light 
curve shape, < 1.1, which are expected to be spectroscopi- 

cally "normal" SNe la (see Figure [141 . 



